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Abstract. The final result is presented of the NA48 experiment, performed at the CERN SPS neu- 
tral kaon beams, on the direct CP violation parameter Re(e'/e), as measured from the decay rates 
of neutral kaons into two pions. The data collected in years 1997-2001 yield the evidence for the 
direct CP violation with Re(e'/e) = (14.7 ± 2.2) x 10~ 4 . Description of experimental method and 
systematics, comparison with the corresponding FNAL result and discussion of some implications 
for theory are given. 
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1. Introduction 

CP violation can occur either via the mixing of CP eigenstates, called indirect CP vio- 
lation, or through the interference of amplitudes with different isospins, called direct CP 
violation. Both processes can be measured in decays of neutral kaons into pions. In one- 
loop approximation, the first process is described by box diagrams and parametrized by 
precisely measured parameter e = (2.28 ± 0.01) x 10 -3 [1], whereas direct CP violation 
is described by penguin diagrams and quantified using parameter e'. Direct CP violation 
is an elusive effect which evaded conclusive measurements in experiments of previous 
generation performed about a decade ago [2]. On the other hand, the standard model of 
electro-weak interactions does not provide with the value of Re(e'/e) precise enough to 
be comparable with experiments. Therefore, it became tempting to perform more accurate 
measurement with high statistics and well controlled systematics in order to conclude on 
the existence of this effect. Two experiments have recently provided results on the mea- 
surement of Re(e'/e): the NA48 at CERN [3] and the KTeV at FNAL [4]. 



2. Experimental method 

Experimentally, it is convenient to measure the double ratio R of frequencies of four si- 
multaneously registered decay modes of neutral kaons into two pions, related to the e' 
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|(A°|K L )/(A°|X S )| 2 

\(n+n-\K L )/(Tr+ir-\K s )\2 
~l-6Re(£'/£) (1) 

thus exploiting cancellations of systematic effects contributing in the same way to all pro- 
cesses. 

In NA48 experiment at CERN, decays of -Kl.s from two collinear beams, produced 
by 400 GeV protons on two targets separated by 120 m (fig. 1 left), are observed in the 
same decay region and registered in detectors downstream (fig. 1 right). Kg decays are 




Figure 1. Schematic view (not to scale) of CERN neutral kaon beamline (left) and the 
NA48 detector (right), located downstream of the decay region. 



distinguished from Kl by a coincidence between the decay time and times of protons 
producing the K$ beam (7r°7r° mode) or by decay vertex position transverse to beams 
(7r + 7r~ mode). In order to minimize difference in acceptance due to different lifetimes, 
Kl decays are weighted as a function of their proper times thus making their spectra 
similar to those of Kg- 

Decays K° —> 7r + 7r~ are measured by a magnetic spectrometer composed of four drift 
chambers with a dipole magnet. A liquid crypton electromagnetic calorimeter (LKr), being 
a hodoscope of 128 x 128 separately read out cells, is used to reconstruct K° — * 7r°7r° 
decays. 

Triggers are designed to reduce the rate of 400 kHz of particles in the spectrometer to 
less than 10 kHz. 

The trigger for K° —> 7r°7r° decays operates on analogue sums of 2 x 8 cells, horizontal 
and vertical, converted into kinematic variables using look-up tables. It requires an electro- 
magnetic energy deposit greater than 50 GeV, a center of energy at calorimeter impact 
smaller than 25 cm, vertex position not exceeding 5 K $ lifetimes and less than 6 peaks in 
each of transverse projections. 

Charged trigger has two levels. At the first level it requires coincidence of three fast 
signals: coincidence in opposite quadrants of the charged trigger hodoscope, chamber hit 
multiplicity and total calorimetric energy deposit above 35 GeV. The second level consists 
of a fast event-building asynchronous processors reconstructing tracks. 

Triggers are accepted if a set of requirements for tracks, vertex and reconstructed mass 
quality are fulfilled. 
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3. Event reconstruction, selections, background rejection and subtracton, acceptance 
corrections and systematic effects 

For K° — > 7r 7T° decays only information from 3-100 GeV photon showers in LKr is used. 
Longitudinal vertex position z vx is reconstructed using shower energies and their impact 
coordinates at LKr. The average resolution on z vx is 55 cm and on the kaon energy is 
0.5%. Invariant masses of 7-pairs are then computed. The 7 pairings are ascribed to 7r°s 
by requiring that the sum of pair masses is closest to double ir° mass and their difference to 
zero. The K° mass is reconstructed from ir°ir° decays with resolution well below 1 MeV. 
Event time is estimated from the most energetic cells in each cluster and is known with 
accuracy of 220 ps. 

The background in neutral mode comes uniquely from Kl — * 37r° decays and is largely 
suppressed by requiring no additional showers within ±3 ns around event time. Remaining 
background is removed by a constrained for \ 2 built from masses of 7 pairs (cf. fig. 2 left). 

The K° — > 7r + 7r~ decays are reconstructed from tracks using hits in the drift chambers, 
the magnetic field map and detectors alignment. Average z vx resolution is 50 cm and that 
for x vx ,y vx is 2 mm. The kaon mass is reconstructed in this mode with accuracy of 2.5 
MeV. Event time is given by charged hodoscope with accuracy of 300 ns. 

Charged background to Ks decays comes from A — > pn~ and is suppressed using 
events with close momenta of the positive and negative particles. For Kl it comes from 
semileptonic decays and is rejected by cut on E/p < 0.8, no hits in muon veto, constraint 
on effective mass and cut on px (fig-2 right). 




Figure 2. The x 2 f° r Kl,s — > 2ty° candidates, where excess in Kl sample for large 
X 2 is due to 3tv° background (left). Comparison of p T for different backgrounds to 
Kl — > tv + iv~ (right) 



For the 7r + 7r~ mode Kl and Ks are distinguished by transverse vertex position. For the 
7r°7r° mode, however, only z vx is reconstructed and one needs time information from the 
tagger (cf. fig. 1 left) and require its coincidence, within ±2 ns, with event time, thus distin- 
guishing Ks from Kl- Systematic errors related to tagging come from tagging inefficiency 
(small probability not to register time coincidence) or mistagging (accidental coincdence 
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between proton and event time). Both effects are estimated using vertex identification from 
charged mode and included to systematic error. 

Since the acceptance of the apparatus depends on z vx and the Kl and Ks lifetimes are 
different, a correction is needed to account for difference in acceptances between the Kl 
and Ks- Distributions of relevant variables for Kl are thus weighted using lifetimes and 
accounting for interference term. The procedure was checked and systematic error was 
calculated using large statistics Monte Carlo (4 x 10 s decays per mode). 

Other systematic effects contributing to the value and errors of R include the collimator 
scatering on the Ks beam, errors in the energy and distance scales from the LKr (inhomo- 
geneities, energy leakage from clusters, non-linearity of energy response, non-Gaussian 
tails in energy response) and accidental effects due to intensity time dependence and inten- 
sity difference between Kl and Ks beams . 



4. The result and discussion 

Overall statistics of reconstructed and accepted events used for determination of R amounts 
to 65 million (cf. the table). 



mode 


events (x 10 b ) 


mode 


events (x 10 b ) 


K L -> 7r U 7T U 


4.7 




— > 7T + 7T 


21.6 


K s -> 7r°7r° 


7.4 


K s 


— > 7T + 7T~ 


31.8 



R stability against cut variations 



Cg< 1 1.5 cm 
Cg<7 cm 



pY<3.0xiar" 
p'f< 1.5x10* 
X 2 < 17.1 
% 2 <9.9 
T < 3.8 t s 
T < 2.9 t s 
tagging window ± 2.5 ns 
tagging window ± 1.5 ns 
accept OX dead time 
accept MBX dead time 
reject DCH-ovfl ± 234 ns 
reject DCH-ovfl ± 350 ns 
no extra clusters in n*7i" 
no Ks/KI intensity weighting 
y radius > 20 cm 
track radius > 20 cm 
mom. asym. < 0.2 
no mom. asym. cut 
ingoing tracks 
outgoing tracks 
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Figure 3. The summary of corrections and systematic errors to 7?, where Rstandard 
corresponds to the final, corrected value of the double ratio, the belt represents system- 
atic error and points with error bars are contributions to the correction on R (left) and 
world data on the Re(e'/ e ) (right). 



Fig. 3 left summarizes all sources of systematic errors and corrections to R. It is seen that 
major systematic effects to R come from the acceptance and charge background. 
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The final value of the CP violation parameter obtained by the NA48 Collaboration is 

Re(e'/e) = (14.7 ± 1.4 ± 0.9 ± 1.5) x 1(T 4 

= (14.7 ±2.2) x (2) 

where the first error is pure statistical from the 2ir samples, the second is systematic from 
the statistics of control sample used for its estimation and the third is the contribution of all 
other systematic errors. Stability of Re(e'/e) as a function of the K° energy was carefully 
checked and no energy dependence was found. Also, NA48 finds good consistency of 
results monitored year by year between 1997 and 2001 . 

World data for Re(e'/e) is presented in fig. 3 right together with the grand average 
from four exeriments [2-4] Re(e'/e) = (16.7 ± 1.6) x 10~ 4 . Good agreement between 
experiments can be judged from the figure. Existence of the direct CP violation in K° 
decays is thus proved, providing with fairly precise number to the theory. It practically 
excludes Wolfenstein's superweak hypothesis [5] which assumed e' = 0. 

As commonly known, expanding e'/e in terms of the matrix elements of the four-quark 
operators Q n (n = 1, . . . , 10), being components of the effective Hamiltonian, one sepa- 
rates it into AI = 1/2 and AI = 3/2 parts [6]. The first part is dominated by gluonic 
penguin process and the second one by the Z° penguin. The latter is proportional to the 
square of the top quark mass and the degree of cancellation between the two, due to their 
destructive interference, is sensitive to it. Since the top quark mass is already known at 
present with accuracy of 5 GeV/c 2 [1] and penguin cancellations are known to be large 
with good precision, larger theoretical uncertainty comes from the long-range part of ma- 
trix elements due to strong interactions which have to be calculated non-perturbatively. 
Therefore current experimental precission for e'/e is not really useful for eventual im- 
provement of the CKM matrix parameter r\ until our ignorance about strong parts obscures 
the whole picture. Efforts towards estimation of QCD factors from lattice calculations still 
do not lead to satisfactory results [7], mainly because of usage of uncontrolled approxima- 
tions, as quenched approximation, in order to make the problem computationally tractable. 
For example, currently obtained values of e'/e are negative and consistent with zero within 
errors, where theoretical errors are larger than experimental [8]. In addition, contribution 
to e'/e from final state interactions is not estimated and it may not be small due to relatively 
large relative momenta of pions in the final state. 
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